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(N| ; Abstract 

The two-body Coulomb Hamiltonian, when calculated in Coulomb-Sturmian basis, has an infinite symmet- 
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ric tridiagonal form, also known as Jacobi matrix form. This Jacobi matrix structure involves a continued 
fraction representation for the inverse of the Green's matrix. The continued fraction can be transformed to a 
ratio of two ^F\ hypergeometric functions. From this result we find an exact analytic formula for the matrix 
elements of the Green's operator of the Coulomb Hamiltonian. 
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I. INTRODUCTION 

Dynamics of a quantum system is fully determined by its Green's operator. The knowledge 
of Green's operator is equivalent to the complete understanding of the system. From the Green's 
operator we can extract the complete spectrum and the wave functions. 

However, in most of the cases, it is not possible to calculate the exact Green's operartor. For- 
tunately, it is often sufficient to determine the Green's operator of an asymptotic Hamiltonian, 
because the rest can be considered as a perturbation and can be approximated by finite matrices. 
This way, having an analytic representation of the Green's operator of some asymptotic Hamilto- 
nian, one can build a powerful quantum mechanical approximation method. 

This is the central idea behind an approximation scheme that is used to solve few-body prob- 
lems. In particular, we have applied this scheme with great success to determine the solution 
three-body Faddeev equations with Coulomb interactions (see eg. 1 1] and references therein). The 
success of the calculation depended very much on our ability to evaluate the matrix elements of 
the Coulomb Green's operator between Coulomb-Sturmian basis states. The choice of Coulomb- 
Sturmian basis is essential. In that basis, the two-body Coulomb Hamiltonian has Jacobi-matrix 
( J-matrix) form. As a result, the matrix elements of the Green's operator satisfy three-term recur- 
sion relations that greatly simplify the calculation of the matrix elements of the Green's operator 
of the Coulomb Hamiltonian. In Ref. B2D we have described two independent ways to calculate 
the matrix elements of the Coulomb Green's operator. They are suitable for different regions of 
the complex energy plane. In the first method, the J-matrix was used as a three-term recursion 
relation. The seed, which is related to 2 F\ hypergeometric function, was derived from indepen- 
dent considerations. In the second method, the explicit inversion of the infinite J-matrix resulted 
in a continued fraction. We also note another method yfl for calculating Green's matrices which 
is based on the solution of scattering problems on L 2 basis J4|]. 

In this article we unify the two approaches presented in Ref. [2]. We start from the continued 
fraction representation for the ratio of two consecutive matrix elements of Coulomb Green's oper- 
ator. A closer inspection reveals that this continued fraction corresponds to the continued fraction 
of the ratio of two 2 -Fi hypergeometric functions. From this we obtain an exact analytic formula 
for the Coulomb Greens matrix. 

In Section II we derive an expression for the inverse of the N x N Green's matrix. A factor 
Cn+i in the inverse of the truncated Green's matrix is expressed as a continued fraction. In section 



Ill, we apply our formulae to a Coulomb problem. We show that Cn is exactly computable. We 
observe that the continued fraction Cn is identical to a certain T-fraction that corresponds to a 
ratio of two hypergeometric functions 2 F\. Establishing the equivalence amounts to solving a 
system of five equations and choosing proper region of convergence of the representation. Finally, 
we establish formulae that are used to compute all elements of the Green's matrix. 

H. CONTINUED FRACTION REPRESENTATION OF THE GREEN'S OPERATOR 



Formally, the Green's operator G is defined by the equation, 



G{z){z-H) = {z-H)G{z) = 1 



(1) 



where z is a complex number and H is the Hamiltonian. 

Suppose that the operator z — H, evaluated in some discrete Hilbert-space basis {\i)}, has an 
infinite symmetric tri-diagonal, i.e. Jacobi-matrix, structure, 
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In such a basis, Eq. © becomes 
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The knowledge of the N x N upper left corner of the full Green's matrix is sufficient to deter- 
mine the physical quantities. Let us denote the corresponding N x N upper left corner matrices 
by j( N \ G^ and 1^, respectively. If we multiply the iVxoo part of J with the ooxiV part of 
G we get the N x N unit matrix 1^). The sum, due to the tridiagonality of J, is reduced to three 
terms 

Jn,n— l"n— X,m T" <-'n,n"n,m ~r , -'n,n+l"n+l,m On,m ) Vv 



where n — 0, 1, ..N and m — 0, 1, ..N. lfn<N, only terms from G^ are appearing in the sum. 
For then = N case, we have: 

JN,N-lGN-l,m + JN,NGN,m + JN,N+lGN+l,m = ^7V,m • (5) 

The Gn+i,™ elements are outside the G^ N > matrix. We can eliminate them formally by writing 

JN,N~\GN-l,m + (Jn,N + JN,N+lG N+l,m/ G N,m)G N,m = &N,m • (6) 

This formal elimination of the elements outside of G^ amounts to modifying single element, 

J N , N Of JW. 

We can calculate the ratio G N+ i^ m /G N ^ m from another, independent relation: 

JN+l,NGN,m + JN+l,N+\GN+l,m + Jn+1,N+2G N+2,m = . (7) 

By rearranging, we get 

1 Gn+i,' 

• (8) 



Jn+i,k G N>m J / 1 Gjv+2)m 

Jn+1,N+1 — Jn+1,N+2 ( j 7; 1 Jn+2,N+1 

JN+2,N+l ( ~*N+l,m 



We introduce a simplifying notation 



1 Gjv+l, m 

CW+i = — = ^ • (9) 

JN+1,N (JN,m 



Then, equation © becomes: 



^ + l = - T = p; = , (10) 



or 



1 

G N+1 = Jn+i,n+i ~ Jn+i,n+2 n -x Jn+2,n+i ■ (11) 

A repeated application of this relation results in a continued fraction. Taking into account that the 
J-matrix is symmetric, we get 

r -l . T J N+l,N+2 



J N+1 ~ J N+1,N+1 

Jn+2,N+2 



T 2 

J N+2,N+3 



Jn+3, 



T 2 



N+3 



This continued fraction does not depend on the index m, the correction term to Jn,n is the same 
for all m's. Consequently, we can write Eq. © in the form 



J, 



ij — Si j NSj,N-JN,N+l^'N+l)G ( — l ( '. 



The modified Jacobi matrix, J, 



(TV) 



i -J 



^i,N^j,NJN n+i^n+i, is the inverse matrix of G^ 



G - (Ju ~ ^i,nSj,nJn,n+iGn+i) 



(13) 



(14) 



III. 7J-DIMENSIONAL COULOMB PROBLEM 



The Hamiltonian of the D-dimensional Coulomb problem, with D > 2, is given by 

h 2 f d 2 L(L + 1)\ Z 



(15) 



2m \dr 2 r 2 J r 

where L = 1 + ^=2-. In the calculation below we set h = m = 1. The Coulomb-Sturmian functions 



are defined by 



ipn{r) 



r(n + l) 



1/2 



.—bar 



{2b s r 



.L+l 7-2L+1 



(26 s r) 



(16) 



T(n + 2L + 2)_ 

where n = 0, 1, • • • , L" is an associated Laguerre polynomial and bs is a parameter. The Coulomb- 
Sturmian functions form a discrete basis, on which the Coulomb Hamiltonian has a J-matrix form 



fl 



•Jn: 



h 2 -b 2 

" °£-{n + L + l)-Z 



for n = m 



k ^^^{n + l){n + 2L + 2) for m = n + l 



(17) 



r + 6 



46c 



£ \/n(n + 2L + 1) for m = n — 1 



where fe = y 7 ^- 

To calculate the Green's matrix of Coulomb Hamiltonian we need to compute the continued 
fraction C^. According to (fTTT) and (IT2T) . the continued fraction Cjv satisfies the relation 
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'-(N + L + l)-Z 



+ (-)^r^(W + l)(N + 2L + 2)C N+1 
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( fk 2 + b 2 
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(JV + L + 1)-Z 



+ K 

p=i 
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(iV + p)(iV + 2L + l+p) 



\ 



^ 2 -h 2 

1 b -£(N + L + l+p)-Z 



26c 



(18) 



From the theory of continued fractions we know that a particular ratio of two hypergeometric 
functions can be represented as a continued fraction from the class of so called T-fractions o] 

T(a, b; c; y) = (c + (b - a + 1) y)+ K 



P =i \c + p+(b-a + l+p)yj 2 F\(a,b+ l;c + l;y) ' 

(19) 

This T-fraction converges to the ratio of two hypergeometric functions if \y\ < 1. Also, for 

y = — 1, the representation is convergent if \$s(c — a + b)\ < |3?(c + a — b — 1)|. 

By comparison, we see that two continued fractions, Eq. (fT8l) and Eq. (fT9l) have identical 
structure. Both fractions have numerators that are quadratic in index p and denominators that are 
linear in index p. Lorentzen and Waadeland [5] have shown that all T-fractions of this type are 
convergent and expressible as a ratio of two hypergeometric functions. 

Hence, we can write, 

tv h \ 2F 1 (a,b;c;y) x 

T(a, b; c; y) = c — — — r = dC N , (20) 

where d is an overall scale parameter. Parameters, a, b, c, y and d are determined from the set 
of five equations subject to a convergence condition \y\ < 1. The five equations that determine 
parameters are: 

*(^f)\ (2D 



.2 _j_ 1.2 ■> :] 



y{b + c-a)= d 2 (^^\ (2iV + 2L + l), (23) 



k 2 + bV ' 



yb(c-a)= rf2|^_r_^j N(N + 2L + 1), (24) 

c+(b-a + l)y= d[-Z+(K^\(N + L + l)\ . (25) 

From the first two equations, together with the convergence condition \y\ < 1, we get 

d = - Ab s (26) 

{bs-ikf 

( b s + ik \ 2 
V= [b^k) ' (27) 

There are two sets of solutions for the remaining parameters, (we also introduce the conventional 



(28) 



parameter 7 = Z/k) 

a = —L + vy 

b = N < 

c = iV + L + l + i 7 

In fact, the two solutions yield the same expression for CV This is easy to show using the identity 
satisfied by the hypergeometric function: 



a = L + 1 + i 7 

b = N + 2L + 1 

c = N + L + l + i>y 



2 F 1 (a, b; c; y) = 2 F X (c - 6, c - a;c;y) x 2 Fi(a, a + b-c;a;y) 
In our case this identity reads, 



2 F 1 (-L + i~f,N;N + L + l + i T ,y) 



2 F 1 {L + l + ij,N + 2L + l;N + L + l + i r , y) 



(i-yf L+1) 

Finally, we have the following closed form expression for the continued fraction 



(29) 



(30) 
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C 



N 



(b S - ik) 



-^ 2^1 



-L + ij,N+l;N + L + 2 + ir, 



bs + ik 
b s -ik 



iV + L + l + i 7 / fbg + ik 

2 F 1 [-L + ijiNiN + L + l + ir 1 



(3D 



b s - ik 

The complete discrete spectrum of the system follows from any matrix element of the Green's 
operator and they are constructed by recursion relations, see equations ® and (ITTT) . In particular, 

Cn+i 
For example, Go,o = C , 
46 S 



— CnJn^-iGn-i,™, 
Jn,n 1 



T 2 

J N,N+1 



J 2 

J N,N+1 






(32) 



G 



1 



0,0 



jFi -L + i 7 ,l;L + 2 + i 7 ; 



bs + ik 
bs - ik 



(33) 



(65 - ik) 2 L + 1 + i 7 

The 2-^1 hypergeometric function has branch-cut singularity on the [1, 00) interval. This cut is 
mapped in Eqs. (l3ll and (l33t to a branch cut in the 2-plane along the positive real axis. In the case 
of the attractive Coulomb potential Z < 0, and i 7 = iZ/k can take negative values. Then, there is 
a first order pole on the negative real axis at L + 1 + i 7 = 0. Additional poles are determined from 
the hypergeometric function by the condition, L + 2 + i 7 = 0, — 1, —2, . . .. All together, there is a 
set of infinitely many discrete first order poles given by the formula, 

* 1 Z 2 



E 



- - 

2 Vi 7 



n r 



2 (n r + L + l) 2 ' 
This is a complete discrete spectrum of the attractive Coulomb system 



0,1,2,... . 



(34) 
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IV. SUMMARY 

In this paper we derived a closed formula for the Coulomb-Sturmian matrix elements of the 
Coulomb Green's operator. In particular, we expressed the inverse of the N x N Coulomb Green's 
matrix in terms of 2 F% hypergeometric functions. The ratio of these two 2 F\ hypergeometric func- 
tion can easily be evaluated in numerical computations by another Gauss-type continued fraction 

a. 
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